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Surface  Correlation  Relationships 
in  Rough  Surface  Scattering 


1.  INTRODUCTION 

There  are  a  large  number  of  approaches  to  the  calculation  of  the  electromag- 

1  _  6 

netic  scattering  from  rough  surfaces.  The  importance  of  these  techniques  rests 

in  their  relevance  to  the  determination  of  the  effects  of  the  real  environment  on 

radar  and  communications  systems.  The  analyses  contain  certain  assumptions 

about  the  nature  of  the  rough  surface  in  relation  to  the  electromagnetic  phenomena. 

Particular  emphasis  has  been  placed  on  the  characterization  of  the  surface  in  terms 

i)f  the  statistical  distribution  of  the  heights  and  their  degree  of  correlation.  These 

features  are  then  related  to  a  normalized  radar  cross  section  for  the  terrain  a- 

7  8 

through  an  electromagnetic  analysis.  ’  We  are  particulars  concerned  with 
questions  involving  the  application  of  physical  optics  principles  to  the  electromag¬ 
netic  analysis.  This  is  important  since  this  approach  is  extensively  vised  both  as 

9 

a  complete  solution  and  as  part  of  composite  models.  “ 

Since  the  surface  height  statistics  and  the  correlation  that  describes  the  height 
relationships  for  pairs  of  points  on  the  surface  are  essential  elements  of  the 
scattering  analysis,  it  is  important  to  understand  how  these  factors  interact  in  the 
analytical  models.  In  this  report,  we  will  describe  several  areas  where  the 

(Received  for  publication  6  March  1986) 

(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 

See  References,  page  27.) 
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statistics  and  correlation  have  important  roles  in  determining  the  scattering  from 
a  random  rough  surface. 

In  the  first  part  of  this  report  we  will  examine  statistical  descriptors  of  a 
surface  and  establish  the  relationship  between  them  and  the  conditions  for  applying 
a  physical  optics  model  to  the  scattering  from  that  surface.  The  second  part  of  the 
report  addresses  the  question  of  how  sensitive  the  distribution  of  scattered  power 
is  to  changes  in  the  form  of  the  surface  height  statistics  and  correlation  function. 
Specifically,  we  will  examine  the  effect  of  three  different  correlation  functions  on 
the  normalized  cross  section  of  a  rough  surface.  For  the  large  scale  physical 
optics  cross  section,  both  Gaussian  and  exponential  surface  height  distributions 
are  considered.  The  small  scale  perturbation  cross  section  is  also  studied  as  a 
function  of  elevation  scattering  angle  for  different  angles  of  incidence  and  different 
levels  of  roughness. 


2.  THE  CONDITIONS  FOR  PHYSICAL  OPTICS 

Over  the  years,  the  arguments  to  support  the  validity  of  physical  optics  have 

centered  around  the  use  of  the  Kirchhoff  integral  representation  for  the  scattered 

em  wave,  where  the  boundary  conditions  on  the  surface  have  been  satisfied  through 

the  use  of  the  Fresnel  plane  wave  reflection  coefficients.  The  boundary  conditions 

are  met  by  specifying  the  total  field  on  the  rough  surface  as  the  sum  of  an  incident 

field  and  a  scattered  field.  The  scattered  field  is  expressed  as  the  product  of  the 

incident  field  and  the  Fresnel  plane  wave  reflection  coefficient.  The  use  of  the 

Fresnel  plane  wave  reflection  coefficients  is  justified  if  the  local  radii  of  curvature 

on  the  rough  surface  are  large  compared  to  a  wavelength  (R^  >  >  A  ).  We  wish  to 

relate  this  constraint  to  statistical  parameters  of  the  surface. 

Utaby  et  al1^  give  two  conditions  which  must  be  satisfied  for  the  Kirchhoff 

2 

approximation  to  be  valid.  These  are  kT  >  6  and  T  /ct  >  2.76  A  where  k  =  2it/A  , 

A  *  em  wavelength.  T  =  surface  correlation  length  and  o  =  standard  deviation  in 
surface  height.  In  this  report,  we  will  demonstrate  that  the  single  condition 

(T  >  N  A  )  is  a  sufficient  condition  to  imply  that  the  radius  of  curvature,  R  ,  must 

c 

be  large  compared  to  a  wavelength  (Rc  >  >  A  ),  for  surfaces  with  a  Gaussian 
distribution  in  heights  and  either  a  power  law  correlation  function  or  a  Gaussian 
correlation  function.  We  will  also  show  that,  except  for  the  case  of  a  surface  with 
small  slopes,  this  condition  is  a  necessary  condition  as  well.  This  means  that,  for 
small  surface  slope  conditions  physical  optics  can  apply  even  when  the  correlation 


10.  Ulaby,  F.  T.,  Moore,  R.K.  ,  and  Fung,  A.  K.  (1982)  Microwave  Remote  Sensin 
Vol.  II,  Addison-Wesley,  Massachusetts. 


length  is  not  large  compared  to  the  wavelength.  However,  if  one  should  want  to 
apply  physical  optics  to  surfaces  with  larger  slopes,  then  the  surfaces  are  limited 
only  to  those  where  T  >  >  A  . 


3.  CORRELATION  LENGTH  AND  RADIUS  OF  CURVATURE 


In  the  previous  section,  we  indicated  that  the  relationships  between  T  and  A 
and  and  A  depended  on  the  magnitude  of  the  surface  slopes  of  the  scattering 
surface.  The  first  point  in  developing  the  specifics  of  these  relationships  is  to 
introduce  the  surface  height  correlation  functions  for  which  the  analysis  has  been 
carried  out.  It  should  be  noted  that  the  surface  is  assumed  to  have  a  Gaussian 
height  distribution  function  and  is  homogeneous  and  isotropic. 

The  surface  is  considered  to  have  either  of  two  types  of  correlation  function. 
The  first  is  the  Gaussian  correlation: 


R(t)  *  a2  exp(-r2/T2) 


where 


(1) 


a  =  standard  deviation  in  surface  heights 

and 

t  =  distance  between  two  points  on  the  surface. 

The  si  rond  form  of  correlation  is  a  power  law; 

■>  ■>  ■> 

R(t)  =  [  1  +  r“/T  ]  P 


where  k  =  a  positive  number.  Cosgriff  et  al'^  used  a  power  law  correlation 

function  with  k  =3/2. 

P 

.Next,  we  introduce  the  expression  for  the  radius  of  curvature  and  relate  it  to 
the  correlation  function  of  the  surface.  The  results  are  obtained  for  each  of  the 
two  correlation  function  types. 


11. 


Cosgriff.  R.  ,  Peake,  W,  ,  and  Taylor,  R.  (  1960)  Terrain  Scattering  Properti 
for  Sensor  System  Design  (Terrain  Handbook  11)~  Ohio  State  L'niversitv, 


PH 


The  average  radius  of  curvature,  R  of  a  Gaussian  surface  described  by  the 
equation  z  =  £  (x,  yi  is  given  by 


1 


<  K  > 


[(1  4-  (z'))2]3/2 


(3 


where  K  is  the  curvature.  The  average  of  (z')  is  given  by 

2, 


2. 


<( z 1 )  >  = 


3£  3£ 

x  Tx 


a  r(t  =  o) 

1 — 


a  t 


■m2 


Also,  the  average  of  (z  )  is  given  by 


<(z  ")2>  = 


a  R(r  =  0) 
1 — 
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For  a  Gaussian  correlation  function, 


<(z')> 


(5 


and 


<(z  > 


12  a 
T4  " 


For  a  powi'r  law  correlation  function, 


iz1  r 


2  k  ct“ 


i  t 


an  I 


<  i  r. 1  ’ 


r  • 


2  2 
12k  ct 


T 


12k  a2 
P 


9 

12  cr“ 


(k  +  1 )k  . 
P  I1 


These  results  relate  the  properties  of  the  slopes  and  slope  derivatives  to 

derivatives  of  the  surface  height  correlation  function.  This  is  consistent  with  the 

theorem  from  random  processes  that  states  that  the  distribution  of  the  derivative 

2 

of  a  normal  process  with  zero  mean  and  variance,  a  ,  is  again  normal  and  has 


a  variance 


(7) 


a 2  =  -a2  p"  (0)  . 
s 

We  next  use  these  results  to  establish  our  relationships.  As  discussed  before, 
Uiaby  et  al^  arrived  at  similar  conditions  under  the  restriction  of  small  slopes. 
We  consider  the  entire  range  of  surface  slopes,  though,  and  hence  can  not  use 
their  simplified  expressions.  As  a  result,  the  analysis  is  more  complicated. 

3.1  General  Solution 

We  have  examined  the  statistical  relations  describing  the  variances  of  the 
surface  height  slopes  and  the  slope  derivatives.  We  now  extend  that  concept  to 
the  covariance  relation  between  these  two  quantities.  As  before,  we  have 

/  2  2 

r  =v/ux)  +  (Ay)  (8t 


and 


>=  3  m+nR(T  > 


3  x 


m+n 


We  now  evaluate 


>  for  the  cases  of  Gaussian  and  power  law  correla¬ 


tion  functions  for  the  surface  heights.  For  the  Gaussian  case,  we  have 


3 

1211(7  )  (x  -x  )  811(7  )  (x-  x  )3 

_p3  3  U(T)  1 

=  /-  n 

O  0 

a  x 

s 

ii 

o 

t 

T4  T6 

=  0  . 


7=0 


(0) 


-Similar  ly,  for  the  power  law  case 


(-  1) 


3  3  "11(7  ) 


3  x 


T 


=  (-1) 


7=0 


12k  (k  -*■  1)11(7  )(x  -  x  )  8k  (k  +  1)11(7  )(x-  x  > 

_ P  ,P  _ „  ,,  °  .  P  £  _ ,,  2 

,_.2.2  * 


T  [1  -  (7  T) 


T  1  1-17  T)  ]' 


!ik  (k  -*  1  )2  11(7  )( x  -  x  )  * 
p  p  o 

xi — 


T[l  +  I7  T)- 


=  0  . 


(10) 


Thus,  for  a  Gaussian  surface  (the  heights  have  a  Gaussian  distribution'  with  either 
a  Gaussian  or  power  law  correlation  among  the  heights,  the  slope  and  slope  deriva¬ 
tives  at  any  point  on  the  surface  are  statistically  uncorrelated.  Further,  the  pre¬ 
viously  cited  theorem  fr0m  random  processes  can  be  applied  to  get  the  result  for 
the  probabilitv  densitv  function  describing  the  slope  at  a  point: 


-1/2 

p(z')  =  |2  jr  <(z')”>]  exp 


Also,  for  the  slope  derivative,  we  have 


J2*  <(z")2>1/2 


2<(z")  > 


Then,  since  z'  and  z"  are  uncorrelated  we  can  write  a  joint  density  function 
p(z',  z")  as 

p(z ' ,  z">  =  |4ir2<(z')2x(z")2>]1/2  exp  -1/2  (  i?-'-*!  +  (z \ 

\<(z')  >  <(z ")  >/ 


Then,  from  Eq.  (3)  we  can  write  the  equation  for  the  expected  value  of  the 
curvature  as: 


<|K|>  =  /  /  I  z"|  Ml  +  <z-)‘ 

-oo  -  00  1  — 


p(z',  z")dz'dz"  . 


Integrating  over  z"  then  yields 

I  I  /  I  9  n  °°  2  ”*^/2  ,  ,,2 

<|k|>  =  -72/77  V<z"Z>  — ■  —  /  [1  +(z')  ]  exp  -^5-i 

•7 2  jt  v  <(z')‘i >  0  2<(z') 


After  some  algebraic  manipulation,  we  obtain 


<  I  K  |  >  =  (V  <(z"l2->/<(z')2>7r)  U  ( 1  /  2 ,  0,  0.  5<  (z')2>  *  1 ) 


where  U(a,  b,  x)  is  the  confluent  hypergeometric  function  of  the  second  kind  with 
x  =  0.5<(z')2>  1  (see  Reference  12). 


12.  Abramowitz,  M. ,  and  Stegun,  I.  A.  (Ed.  )( 1964)  Handbook  of  Mathematical 
Functions  With  Formulas,  Graphs  and  Mathematical  Tables,  N&S  ADDlied 


athematics  series 


1  >IL 


3.2  Different  Slope  Regimes 


Now,  having  obtained  this  analytic  expression  for  <  !  K.  i  >  (or  K  ),  we  want 


examine  the  relationships  between  Rc  and  a  and  T  for  various  slope  regimes. 


As  a  first  step  we  examine  the  small  slope  regime  using  the  general  solution. 
An  asymptotic  expansion  of  U(a,  b,  x)  in  terms  of  the  inverse  of  the  slope  leads  to 
a  zero  order  solution 


<|k|>=  (  '/<(z")2>/n  I  (n/<(z')2>)  (  1  /  (2<  (z' )2  > 


,,  2 

2<(z  )  >/; 


The  relations  for  <(z")  >  for  the  two  correlation  functions  then  give: 


<  !  K  |  >  =  2.  76  a/ T 


(Gaussian) 


<  I K  I  >  *  2.  76  o  ^kp( kp  +  1)/T 


(Power  Law)  . 


These  results  are  equivalent  to  the  ones  obtained  by  Ulaby  et  al  for  small  slope. 


Then,  from  o/T  «  1.  we  have  Rc  >>  0.  4T  for  Gaussian  cases  and 

R  »  0.4T/"7k  (k  +  1)  for  power  law  correlation.  These  results  yield  the  suffi- 
c  p  p 

cient  condition  that  T  »  A  implies  R  »  A  for  Gaussian  correlation  and  for  any 

c  11 

power  law  case  where  kp  ~  1  (including  the  3/2  law  used  by  Cosgriff  et  al.  ).  One 


particular  point  should  be  made  here.  The  case  where  T  =  A  also  leads  to  the 
results  R  »  A  for  the  Gaussian  case  (and  a  similar  argument  would  apply  for  some 
power  law  cases).  Thus,  for  the  small  slope  regime  we  have  the  overall  result 
that  T  >>  A  is  indeed  a  sufficient  condition  for  physical  optics  models  to  apply  but 
it  is  not  a  necessary  condition.  As  long  as  the  surface  slopes  are  small,  even 
surfaces  with  only  a  small  degree  of  correlation  (T  =  A)  can  satisfy  this  condition. 

For  small  slopes,  the  asymptotic  expansion  of  U(a,  b.  x)  to  first  order  gives 


the  results: 


<  |  K  |  >  =  |2.  760/T2l  [1  -  3/2«7/T)2] 


for  Gaussian  correlation,  and 


<  i K |  >  =  (2.  76  a/T2l  4 k  (1  +  k  )  ( 1  -  3/ 2  k  (a/T)2] 

p  P  P 


for  Power  Law  correlation. 


•  .  v./tvL-'.v v -v-  '• 
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We  have  already  showed  the  connection  between  T  >>  A  and  R^  »  A  for  the 
zeroth  order  solution.  For  the  Gaussian  first  order  result  we  have 

R  a  (T2/2.  76  a)  (1  +  3/2(ct/T)2) 
c 

or 

R  >  T2/(2.  76  0)  =S>  R  »  T 
c  c 

and  hence 


T  »  A  R  »  A  for  CT/T  «  1  . 

c 


(21) 


The  analysis  is  similar  for  the  power  law  case.  Thus,  the  general  solution  gives 

T  >>  A  as  a  sufficient  condition  for  R  »  A  for  both  correlations  when  o/T  <<  1. 

c 

As  in  the  zero  order  result,  T  a  A  also  gives  Rc  >>  A  so  the  condition  is  still  not 
a  necessary  one. 

The  next  solution  for  the  general  result  is  that  in  the  large  slope  regime;  the 
intermediate  slope  region  is  deferred  to  the  end  since  the  simplifying  solutions 
do  not  apply  there. 

2  - 1 

For  large  slope  conditions  (CT/T  >>  1),  (2  <z'  >)  -*  0  and 


<  I  K  I  >  =  (  2  /  7T ) 


Uz',2>l<z'2>  +  0[(1/2<z'2»  In  (1  /  2<z'2»] 


We  now  introduce  the  results  for  the  two  correlations  and  get 


R  ~  0.6T  (Gaussian)  and  R  ~  0. 6T/Vk  (Power  Law)  .  (23) 

c  c  p 

As  a  result  we  have  that  T  »  A  is  both  sufficient  and  necessary  for  R  >>  A  for  the 

Gaussian  case  where  a/T  >>  1.  In  addition,  as  long  as  the  restriction  k  =  1  holds, 

P 

this  is  also  true  for  power  law  correlation  functions. 

Finally,  we  address  the  most  complicated  case  where  we  have  to  evaluate  the 
function  explicitly  for  various  intermediate  slope  cases.  These  generally  corres¬ 
pond  to  intermediate  values  of  the  argument  X  but  there  is  some  overlap  into  the 

asymptotic  regimes.  Procedures  for  evaluating  the  function  U(  1  / 2,  0,  X)  are 

12 

outlined  in  Abramowitz.  To  see  how  R  varies  in  the  intermediate  X  range  we 


8 


consider  several  cases  as  10.  0  =  X  -  0.1.  Over  this  range,  for  the  Gaussian 

case.  1.5  2  (o/T)  ^  0.  16,  and  the  evaluation  of  U(l/2,  0,  X)  leads  to  the  result 

0.  5T  5  i  2.5T.  To  complete  the  correspondence  between  intermediate  ranges 

of  slopes  and  X-values  we  must  include  the  small  X  solution  for  0.  1  =  X  2  0.  0025. 

Here  Rc  =  0.  6T.  Similarly,  the  large  X  solution  for  1.0SO/TS  10  results  in 

2.  3T  S  Rc  =  3.6T.  The  remaining  intermediate  range  results  are  summarized  in 

Table  1.  Thus,  over  the  entire  range  of  intermediate  o/T  values,  R^  is  of  the 

same  order  as  T  and  hence  in  that  range  T  »  A  ■*  R  »  A  and  R  A  ■*  T  >>  A. 

c  c 

Thus  T  »  A  is  both  necessary  and  sufficient  for  R^  »  A  when  we  have  intermediate 
o/T  values  and  Gaussian  correlation. 


Table  1.  Results  for  Rc  in  the 
Ilntermediate  Range  of  o/T  Values 


<  o/T  <|k|>  r 


For  the  power  law  case  the  tabular  results  are  similar  with  o/T  including  a 


vk^-term  and  R^  having  a  +  1-term.  Thus,  as  long  as  k^  =  l,  the  necessarv 
and  sufficient  relationship  applies  to  this  correlation  case  as  well. 

At  this  point  we  have  established  that  there  is  a  direct  relationship  between 
the  correlation  length  of  a  surface  and  the  radius  of  curvature  as  far  as  the  mstifi 
cation  for  the  use  of  a  physical  optics  model  is  concerned.  We  have  shown  that, 
except  for  small  slopes,  the  class  of  sui'faces  for  which  R  >>  A  is  restricted  to 
those  for  which  T  >>  A  . 

One  additional  comment  is  worth  mentioning.  It  is  well  known  that  as  the 
average  surface  slopes  become  large  (ff/Tg  1),  the  physical  optics  current  does 
not  truly  describe  the  behavior  of  the  surface  currents.  Shadowing  becomes 
important  under  those  conditions.  The  traditional  approach  to  the  current  is  to 
modify  the  cross  section  obtained  from  physical  optics  by  a  shadowing  correction 


factor  S,  which  describes  how  much  of  the  surface  is  unlit,  that  is,  where  J  -»  0  in 
the  shadow  regions.  12  There  has  been  discussion  as  to  the  validity  of  this  approach. 

Having  shown  that  the  surface  height  correlation  is  an  essential  criterion  for 
establishing  the  validity  of  applying  physical  optics  principles  to  the  analysis  of 
rough  surface  scattering,  we  now  turn  to  the  second  theme  of  the  report.  We 
examine  how  use  of  particular  surface  statistics  and  correlation  functions  affect 

14 

the  actual  determination  of  the  scattering  for  a  range  of  surface  roughness  levels. 


4.  CORRELATION  FUNCTION  AND  SCATTERING 

In  this  section  we  are  concerned  with  the  changes  that  occur  in  the  distribution 
of  the  scattering  from  a  rough  surface  as  a  function  of  roughness  when  alternative 
forms  of  surface  height  statistics  and  correlation  function  are  used  in  the  calcula¬ 
tion.  The  study  is  carried  out  in  terms  of  the  normalized  scattering  cross  section 
of  the  surface,  o’.  The  variation  of  o’  with  angle,  roughness  and  correlation 
function  is  examined  for  different  surface  height  distributions.  The  primarv 
concern  is  the  regime  where  T  »  A  and  physical  optics  models  can  be  used.  For 
completeness,  the  results  when  the  surface  heights  and  correlation  lengths  are 
appropriate  to  perturbation  method  analysis  are  also  considered. 

Three  basic  forms  for  the  correlation  function  are  assumed  in  this  section; 

(1)  The  Gaussian  correlation  function,  given  by  Eq.  (1), 

(2)  The  Bessel  function  correlation  function 

R(T  )  =  a2[l  +  1/8(T  / T)2](t  /T)Kj(t  IT)  -  C2(T  /T)2Kq(T  IT),  (24) 

(3)  The  power  law  correlation  function,  given  by  Eq.  i2). 

Because  of  its  tractability  in  analytical  expressions  for  o°,  the  Gaussian  form 

has  been  widelv  used  bv  many  authors  over  the  years  (Beckmann  and  Spizzichino,  1 
2  j  o 

Ruck  et  al,  Ulabv  et  al  ).  However,  not  all  terrain  nor  sea  states  can  be  repre¬ 
sented  accurately  by  a  Gaussian  correlation  function,  because  the  Gaussian 
correlation  function  has  a  finite  dc  component  in  the  power  spectral  domain  (at 
k  0  where  k  =  wavenumber).  The  Bessel  function  R(T  ',  form  has  been  used 


13.  Sancer,  M.I.  (1969'  Shadow-corrected  electromagnetic  scattering  from  a 
randomly  rough  surface,  IEEE  Trans.  Antennas  and  Propag.  , 

AP-  17 :577-585. 


14.  Brown,  G.S.  (1984)  The  validity  of  shadowing  and  corrections  in  rough 
surface  scattering.  Radio  Science  14,  6:1461-1468. 

-  ^ 


15.  Miller,  L.S.,  Brown,  G.S.  ,  and  Hayne,  G.  S.  (10721  Analysis  of  Satellite 
Altimeter  Signal  Characteristics  and  Investigation  of  Sea -Truth  Data 


where 


j  =  [  i»  1  <2ti>  T  J(V  T)[(l+  2/3  2?(1  -c(t)))‘3/2  -  (l  +  2/3rV3/2]7d 

u  lx/  0  0  xy 


[(1+2/3  22)'3/2  -  (1  +  1/3  £2)"3]  6(vSPEC)  6(vSPEC) 
x  +  x  y 


Here.  A  =  unit  area  of  rough  surface.  The  term  given  by  Eq.  (32)  represents 

the  diffuse  power  scattered  into  all  directions,  whereas  the  J,  term  given  by 

Eq.  (33)  represents  an  additional  amount  of  diffuse  power  scattered  into  the 

1 6 

specular  direction  (Brown  ). 

These  expressions  represent  the  general  form  of  the  solution.  In  order  to 
examine  the  effects  of  using  different  types  of  correlation,  we  must  replace  c(T) 
in  Eq.  (31)  and  Eq.  (32)  by  the  three  explicit  forms  that  we  wish  to  study.  These 
are  given  by  Eqs.  (1),  (2),  and  (24).  In  the  next  section  we  will  show  how  the 
successive  correlation  functions  affect  the  roughness  dependence  of  surfaces 
with  Gaussian  height  distributions.  Then,  we  will  consider  the  effect  of  surface 
type  bv  comparison  with  similar  results  for  surfaces  with  exponential  height 
distributions. 

4.2  Correlation  Function  and  Surface  Roughness 

In  this  section  we  present  results  for  the  changes  in  rough  surface  scattering, 
as  surface  roughness  increases.  The  effect  of  using  different  correlation  functions 
is  shown.  The  roughness  is  introduced  in  terms  of  the  Rayleigh  roughness 
parameter,  Z.  The  surfaces  are  chosen  to  satisfy  physical  optics  conditions.  The 
results  are  not  limited  to  any  particular  surface  type  since  the  dependence  on 
dielectric  constant  has  been  excluded  by  making  comparisons  in  terms  of  the 
quantity  .1  defined  in  Eq.  (29).  Results  are  given  for  azimuth  angle  0g  =  0C.  The 
relative  breadth  of  the  scattering  pattern  is  shown  by  the  amount  scattered  for 
successive  values  of  8  g.  The  surface  correlation  length  is  taken  to  be  T  =  13  m. 
The  wavelength  is  X  =  0.  275  m.  The  incident  elevation  angle  is  ft  =  85°.  Although 
the  results  presented  in  this  report  are  all  at  this  incident  angle,  additional  calcu¬ 
lations  were  performed  at  small  and  intermediate  6h,  and  indicate  that  the  results 
are  similar  for  all  cases.  Then  the  results  are  given  for  Z ^  corresponding  to 
<7.  =  3.  16  m;  for  <7„  =  0.  316  m;  and  Z  for  (7-  =  0.  0316  m. 
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4.  2.  1  RESULTS  FOR  GAUSSIAN  SURFACE  HEIGHTS 

Table  2  shows  the  =-  ttering  cross  section  pattern  for  a  Gaussian  surface  for 
large  Rayleigh  parameter,  Ej.  Here,  Cj  is  the  Gaussian  correlation  case,  Cg  is 


the  Bessel  case  and  Cg  is  the  power  law  case.  The  table  shows  a  broad  pattern 


increasing  in  magnitude  at  large  9  for  all  three  forms.  The  Gaussian  and  power 


law  results  are  similar  and  tend  to  exceed  the  Bessel  case  except  at  small  6  . 


Table  3  shows  the  results  for  intermediate  Rayleigh  parameter  Eg.  Here,  the 
overall  pattern  has  narrowed  somewhat  but  the  trend  is  the  same  as  in  Table  2. 


Table  4  shows  the  results  for  small  Rayleigh  parameter,  Eg.  Here  the  contribu¬ 


tions  are  mainly  in  directions  close  to  specular. 


Table  2.  Results  for  Gaussian  Surface  With  Large  Roughness 


m 

& 


m 


Table  3.  Results  for  Gaussian  Surface  With  Intermediate 
Roughness 


6 

s 

JGA(cl) 

JGA(C2) 

JGA|C3> 

0 

<  10‘7 

0.  008 

<  0.  1 

5 

<  10"7 

0.  011 

<  0.  1 

10 

<  10"7 

0.  016 

<  0.  1 

15 

<  10"' 

0.  024 

<  0.  1 

20 

<  10"7 

0.  038 

<  0.  1 

25 

<  10"7 

0.  063 

<  0.  1 

30 

<  10"7 

0.  108 

<  0.  1 

35 

<  10"7 

0.  196 

<  0.  1 

40 

<  10"7 

0.  382 

<  0.  1 

45 

<  10"7 

0.  801 

<  0.  1 

50 

<  10"7 

1.  828 

<  0.  1 

55 

<  10"7 

4.  595 

<  0.  1 

60 

61  X  10"7 

12. 889 

<  0.  1 

65 

0. 008710 

41. 187 

0.  6 

70 

3.  146 

155. 339 

27.  1 

75 

308.  6 

743. 807 

657.  9 

80 

9313. 0 

5,501. 954 

8,  649.  0 

85 

60.  690.  0 

10,464.87 

97, 600.  0 

90 

23,  245. 0 

18,  948.  55 

19,  450.  0 

Combining  these  three  sets  of  results  confirms  that  for  all  three  correlation 
functions  the  patterns  are  similar  to  those  of  Beckmann  and  Spizzichino.  For 
smooth  surfaces,  scattering  tends  to  be  in  the  specular  direction  and  the  scattering 
pattern  broadens  as  the  roughness  level  of  the  surface  increases.  It  should  be 
noted  that  the  specular  point  null  for  Bessel  correlation  in  Table  3  is  present  only 
for  9  .  =  85°.  There  is  no  similar  result  at  other  6  .  values. 


Table  4.  Results  for  Gaussiar  Surface  With  Small  Roughness 


,JGA(cl) 

JGA(C2) 

JGA(c3> 

0 

A 

O 

i 

-4 

<  0.  001 

<  0.  1 

5 

A 

O 

I 

-4 

<  0. 001 

<  0.  1 

10 

<  io'7 

<  0.  001 

<  0.  1 

15 

t- 

1 

O 

V 

<  0.  001 

<  0.  1 

20 

<  io'7 

<  0.  001 

<0.  1 

25 

<  io'7 

<  0.  001 

<  0.  1 

30 

<  io-7 

<  0.  001 

<  0.  1 

35 

r- 

O 

V 

0.  001 

<  0.  1 

40 

<  io'7 

0.  002 

<  0.  1 

45 

t— 

i 

O 

V 

0.  005 

<  0.  1 

50 

<  io-7 

0.010 

<0.  1 

55 

<  io'7 

0.  024 

<  0.  1 

60 

<  io'7 

0.  064 

<  0.  1 

65 

<  10‘7 

0.  214 

<  0.  1 

70 

<  10-7 

0.  968 

0.  1 

75 

0. 000367 

7.464 

0.6 

80 

95.  78 

166. 792 

168.  0 

85 

1841.  0 

5.  000 

3673. 0 

90 

304.  0 

302.  896 

252.  1 

4.2.2  RESULTS  FOR  EXPONENTIAL  SURFACE 
HEIGHTS 

In  Eq.  (31)  we  indicated  that  the  equivalent  surface  term  for  the  exponential 
case  had  two  elements.  The  Jg  component,  defined  in  Eq.  (33),  represents  an 
incoherent  scattered  contribution,  solely  in  the  specular  direction.  As  can  be  seen 
in  Eq.  (33),  the  term  is  independent  of  correlation  function.  The  magnitude  of  Jg 
was  calculated  for  our  three  Z  values,  and  for  all  cases  it  was  orders  of  magnitude 
less  than  the  associated  term  for  all  three  forms  of  the  correlation  function. 
This  result  also  applies  to  other  incident  angles.  As  a  result  of  this  finding  we 
assumed  that  =  and  made  our  comparisons  using  the  results  for  J  . 

Since  we  wish  to  make  additional  comparisons  in  this  section,  the  tables  of 
results  are  more  complicated  than  the  previous  set.  The  first  comparison  is  the 
equivalent  roughness  effect  for  the  three  correlation  functions.  The  second  aspect 
is  comparison  of  the  results  for  the  Gaussian  surface  height  distributions  with 
those  for  the  exponential  surface.  Table  5  shows  the  corresponding  patterns  for 
both  surfaces  for  a  Gaussian  correlation  function  and  all  three  levels  of  roughness. 
Table  6  presents  the  same  cases  for  a  Bessel  function  correlation  and  Table  7 
illustrates  the  results  for  a  power  law  correlation  function. 

Consideration  of  the  three  tables  shows  that  the  broadening  of  the  scattering 
pattern  as  the  surface  roughness  increases  for  the  exponential  height  cases  is 
clearlv  present  for  all  three  correlation  functions.  The  largest  values  occurred 
for  Eg*  for  all  cases.  The  power  law  and  Gaussian  results  tended  to  be  equivalent 
and  slightly  higher  than  the  Bessel  function  correlation. 

The  comparisons  for  both  surfaces  with  a  given  form  of  correlation  also  show 
very  similar  results  for  both  cases.  There  appears  to  be  a  slight  dominance  for 
the  exponential  surface  at  the  highest  roughness  condition  but  overall  the  scattering 
does  not  appear  to  be  very  sensitive  to  either  correlation  function  or  surface 
height  statistics. 


Scattering  Pattern  Comparisons:  Gaussian  Correlation  Function 


G95  1.0  20,  259.  303.  I  3007.  23,  245.  304. 


ble  6.  Scattering  Pattern  Comparisons:  Bessel  Correlation  Funct 


4.3  Perturbation  Regime 


The  preceding  results  and  discussions  are  all  for  the  case  of  a  physical  optics 
model  where  the  heights  are  large  compared  to  the  wavelength.  The  third  topic- 
in  this  report  is  the  extension  of  these  considerations  to  the  small  scale  height 
regime  (heights  small  compared  to  the  wavelength)  where  a  perturbation  method 
scattering  cross  section  applies. 

4.  3.  1  PERTURBATION  METHOD 

In  this  regime  we  have  an  alternative  formulation  for  the  normalized  scattering 

cross  section  cr  .  The  conditions  under  which  this  relation  can  be  applied  are 
s  s 

kff  <  1  (34) 

s 


and 


a  IT 
s  s 


<  1 


where  a  is  the  fine  scale  height  variation,  and  T  represents  a  correlation  length 
for  the  small  scale  heights.  The  solution  for  <*ass  then  is  given  by  Ruck  et  al  ; 


cr° 


ss 


(35) 


with 


4  2  2  2  i  ! 

J  =  (2tt  )  f(4k  /n)  a  cos  0.  cos  0  |  /  c(t)  J  (  v  T)  rd  v  .  (36) 

ss  1  S  i  s  0  o  xv 

where 

a  =  scattering  matrix  elements  for  small  scale  of  roughness 
[see  Ruck  et  al^]. 

Note  that  the  value  of  CT<,SS  is  independent  of  the  surface  height  statistics  and 
depends  only  on  the  form  of  the  correlation  function.  Thus  the  surface  dependent 
results  of  the  physical  optics  conditions  do  not  apply  here,  and  we  are  concerned 
only  with  the  differences  for  the  three  types  of  correlation  function  as  given  in 
Eqs.  (1),  (2),  and  (24).  The  resulting  integrals  for  Jgs[c(TM  can  be  solved 
analytically  and  we  have 


(1)  Gaussian  correlation  (c^) 


Jgs  =  [<4/tt)k4  a2  cos2  9 .  cos2  0g]  [ttT2  exp  {-T2  v2/4  { 


(37) 


(2)  Bessel  function  correlation  (Cg) 


J  =  [4k2  T  a  cos  9.  cos  9  ]2 

SS  1  S  S  1  SJ 


3  v4  T4 
X  s 


- 2 — 2~T 

(i  +  vx  Ti> 


(38) 


(3)  Power  law  correlation  (c^) 


Jss  =  [81<4  al  cos2  9i  COs2  0s^  Ts  exp  I'  !vx^  Tsi  ‘  (39) 

The  variations  of  a°  were  studied  in  the  same  context  as  was  done  for  the 
ss 

physical  optics  cases.  We  eliminated  the  dependence  on  the  scattering  matrix 

elements  and  make  comparisons  for  the  three  forms  of  J  [  c(T )] . 

s  s 

4.3.2  PERTURBATION  RESULTS 

As  before,  we  chose  A  =  0.  275  m  and  9 .  =  85°.  Here  T  =  0.  0444  m  and  we 

l  s 

considered  three  levels  of  roughness.  k<7g  =  0.  914,  0.  5,  and  0.  2  m.  The  results 
are  presented  in  Tables  8,  9  and  10,  respectively.  The  behavior  is  quite  distinct 
from  that  of  the  large  scale  solutions.  All  three  levels  of  roughness  show  similar- 
angular  distributions  of  scattered  power  and  we  will  discuss  the  results  together. 

The  Bessel  function  correlation  showed  somewhat  different  behavior  from  those 
of  the  power  law  and  Gaussian  cases.  In  all  correlation  cases,  the  relative  magni¬ 
tudes  decreased  with  decreasing  roughness.  There  is  no  tendency  to  peak  near  the 
specular  direction  as  was  true  for  large  scale  roughness.  In  contrast,  all  peaked 
at  small  0g  with  the  Bessel  function  results  being  monotonic  and  the  power  law 
results  always  peaking  at  a  slightly  higher  value  of  0  .  The  location  of  the  peak 
was  independent  of  roughness.  The  Gaussian  and  power  law  results  have  about  the 
same  magnitude  (power  law  slightly  higher)  and  the  spread  of  the  scattered  power 
for  both  is  fairly  uniform  up  to  0g  =  60°.  The  Bessel  correlation,  though,  drops 
off  rapidly  with  scattering  angle  0g  and  can  be  an  order  of  magnitude  below  the  peak 

when  9  reaches  60°. 
s 


Table  8.  Small  Scale  Scattering  Angular  Dependence, 

kcr  =  0.  914 
s 


*s 

J3S<=i» 

Jss<c2> 

Jss(c3 

0 

0, 202629E-01 

0.  196216E-01 

0.  1903339-01 

5 

0. 209883E-01 

0.  192294E-01 

0.  206372E-01 

10 

0. 213185E-01 

0.  178451E-01 

0.  220193E-01 

15 

0.  212233E-01 

0.  155230E-01 

0.  230965E-01 

20 

0.  206946E-01 

0.  1251 15E-01 

0.  237859E-01 

25 

0. 197475E-01 

0.  9226 13E-02 

0.  240126E-01 

30 

0. 184199E-01 

0.  613994E-02 

0.  237  177E- 01 

35 

0. 167696E-01 

0.  363386E-02 

0.  228657E-01 

40 

0.  148701E-01 

0.  188213E-02 

0.  214528E-01 

45 

0.  128057E-01 

0.  83747 0E-03 

0.  195123E-01 

50 

0.  106658E-01 

0.  312801E-03 

0.  17 1 184E-01 

55 

0.  854008E-02 

0.  949424E-04 

0.  143856E-0 1 

60 

0.651372E-02 

0.  223534E-04 

0.  1 14645E-0 1 

65 

0.  466422E-02 

0.  376469E-05 

0.  85325  IE-02 

70 

0.  305867E-02 

0.  392762E-06 

0.  578100E-02 

75 

0. 175258E-02 

0.  187  102E-07 

0.  339986E-02 

80 

0. 789067E-03 

0.  174691E-09 

0.  156004E-02 

85 

0.  198783E-03 

0.  140398E-  1 1 

0.  397566E-03 

90 

0. 

0. 

0. 

1.4  (’.<>iii|>< isilo  Surface 

The  small  scale  surface  results  and  the  large  scale  height  results  have  been 

presented  as  separate  cases.  As  a  final  comment  it  should  be  noted  that  actual 

2 

surfaces  are  often  combinations  of  both  types  of  height  conditions.  Ruck  et  al 

have  discussed  how  a  composite,  two  scale  surface  may  be  obtained  from  the  two 

solutions  (as  long  as  the  slopes  at  the  large  scale  roughness  are  small,  a/ T  <  1). 

The  cross  section  of  the  small  scale  roughness  [Eq.  (35)]  is  simply  added  to  that 

of  the  large  scale  roughness  [Eq.  (29)].  If  the  large  scale  surface  slopes  are  not 

small,  then  the  cross  section  for  a  composite,  two  scale  surface  is  more  compli 

17 

cated.  This  is  discussed  by  Brown. 


17.  Brown,  G.S.  (1978)  Backscattering  from  a  Gaussian  distributed  perfectly 
conducting  rough  surface,  IEEE  Trans.  Antennas  and  Propag. 

APS-26: 47  2-482. 
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Table  9.  Small  Scale  Scattering  Angular  Dependence, 

kcr  =0,5 
s 


J  <c.) 
ss  1 

Jss(c2) 

Jss(c3> 

0.  606388E-02 

0.  587 194E-02 

0.  569606E-02 

0.  628094E-02 

0,  575458E-02 

0.  6 17586E-02 

0.  637976E-02 

0.  53404  IE-02 

0. 658952E-02 

0.  635  126E-02 

0.  46454  IE-02 

0.  691 187E-02 

0.  619303E-02 

0.  3744 18E-02 

0.  7  1 1817E-02 

0. 590963E-02 

0.  276100E-02 

0.  7  18604E-02 

0.  551 236E-02 

0. 183744E-02 

0.  709776E-02 

0.  501847E-02 

0.  108747E-02 

0. 684280E-02 

0.  445002E-02 

0.  563243E-03 

0.  64 1996E-02 

0. 383221E-02 

0.  250620E-03 

0.  583926E-02 

0. 3 19183E-02 

0,  936086E-04 

0.  512284E-02 

0. 255569E-02 

0.  284281E-04 

0.  430504E-02 

0.  1 94930E-02 

0.  668948E-05 

0.  343087E-02 

0.  139581E-02 

0.  1 12662E-05 

0.  255343E-02 

0.  9 15335E-03 

0.  1  17538E-06 

0.  173002E-02 

0.  524477E-03 

0.  5599 1 8E- 08 

0.  101744E-02 

0.  236  135E-04 

0,  522780E-10 

0.  466857E-03 

0.  59487 8E-04 

0.  420153E-12 

0.  1 18975E-03 

Table  10.  Small  Scale  Scattering  Angular  Dependence, 


k(7 

s 

=■  0.  2 

9 

s 

.1  (c  ) 
ss  1 

J  <C„) 

ss  2 

Jss(c3> 

0  0.  970216E-03  0.  93951  IE-03  0.  91  1373E-03 

5  0.  100495E-02  0. 920733E-03  0. 988138E-03 

10  0.  102076E-02  0.  854467E-03  0.  105432E-02 

15  0.  101620E-02  0.  743264E-03  0. 110590E-02 

20  0.  990889E-03  0.  599070E-03  0.  1  13890E-02 

25  0.  945540E-03  0.441760E-03  0.  114976E-02 

30  0.  881973E-03  0.  293990E-03  0.  U3564E-02 

35  0. 802953E-03  0.  173995E-03  0.  109485E-02 

40  0.  7 12005E-03  0.90U94E-04  0.  102719E-02 

45  0.  613154E-03  0. 400993E-04  0.  934280E-03 

50  0.  510693E-03  0.  149774E-04  0.  819653E-03 

55  0. 40891  IE-03  0.454849E-05  0.  688806E-03 

60  0.311 888E-03  0.  107032E-05  0. 548939E-03 

65  0. 223330E-03  0. 180259E-06  0. 408549E-03 

70  0.  146454E-03  0.  188259E-07  0. 276804E-03 

75  0. 839 16 1 E -04  0.  895869E-09  0.  162790E-03 

80  0. 377816E-04  0.836447E-11  0. 746971E-04 

85  0. 951799E-05  0.  672242E-13  0.  190360E-04 


5.  CONCLUSIONS 


t 

i 

This  report  has  been  concerned  with  the  relationship  between  surface  height  \ 

statistics  and  correlation  and  the  scattering  of  em  waves  by  the  surface.  The  first  j 

analysis  showed  that  establishing  conditions  on  the  correlation  length  (T  »  A.)  was 
sufficient  to  cause  the  radius  of  curvature  to  be  large  compared  to  the  wavelength 
as  well.  Hence  this  allows  us  to  establish  the  conditions  for  which  physical  optics 
principles  can  be  used  to  analyze  the  scattering  from  the  random  rough  surface. 

The  second  major  theme  of  the  report  examined  the  manner  in  which  the  pattern  of 
scattered  power  as  a  function  of  the  surface  roughness  is  dependent  on  the  surface 
height  statistics  and  the  form  chosen  for  the  correlation  function  of  those  surface 
heights. 

The  first  topic  involved  the  determination  of  limits  for  the  application  of  the 
physical  optics  scattering  cross  section  model.  The  first  restriction  in  the  analysis 
is  that  the  surface  height  distribution  be  Gaussian.  Under  this  constraint  it  was 
shown  that  the  surface  slopes  and  slope  derivatives  are  statistically  uncorrelated. 

Then,  for  surfaces  with  either  a  Gaussian  or  power  law  correlation  function,  we 
conclude  that  for  intermediate  or  large  slopes  the  condition  T  »  A  is  both  neces¬ 
sary  and  sufficient  for  Rc  »  A.  That  is,  if  one  wishes  to  apply  physical  optics  in 
those  regimes,  one  is  restricted  to  the  class  of  surfaces  for  which  T  »  A.  On  the 
other  hand,  for  small  slope  regimes  the  applicability  is  less  restrictive.  There, 

T  »  A  is  still  a  sufficient  condition  for  Rc  »  A  but  it  is  no  longer  necessary. 

Surfaces  whose  correlation  length  is  on  the  order  of  wavelength  (T  *  A  )  can  also 
generate  the  condition  Rc  >>  A  .  A  further  aspect  of  this  is  that  since  either  T  or  A 
can  be  varied  there  is  no  particular  restriction  on  the  value  of  the  Rayleigh 
parameter  in  the  various  slope  regimes. 

In  the  second  topic  addressed  in  this  report,  a  study  is  made  of  three  aspects 
of  the  angular  dependence  of  rough  surface  scattering.  The  results  assess  the 
dependence  on  surface  height  statistics  and  surface  correlation  function,  but  do 
not  include  the  effects  of  surface  dielectric  constant  in  the  scattering  matrix 
component  of  the  scattering  cross  section.  For  physical  optics  conditions,  the 
results  show  that  for  a  Gaussian  surface  height  distribution,  the  scattering  angle 
dependence  for  all  three  correlation  functions  is  similar.  For  a  smooth  surface 
the  contributions  are  mainly  in  the  specular  direction  and  the  pattern  becomes 
broader  as  surface  roughness  increases.  The  Gaussian  and  power  law  correlation 
results  are  about  the  same,  while  the  Bessel  correlation  results  in  more  variability. 

The  magnitude  is  greatest  for  intermediate  roughness.  Next,  the  dependence  of 
surface  statistics  was  examined.  For  the  case  of  exponential  surface  heights, 
there  is  an  additional  diffuse  contribution  in  the  specular  direction,  independent  of 
correlation  function;  for  the  cases  considered,  this  term  is  always  less  than  the 
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general  diffuse  scattering  term.  The  resulting  distributions  repeated  the  roughness 

dependence  of  the  Gaussian  surface.  Also,  the  Gaussian  and  power  law  results 

are  similar  and  the  Bessel  correlation  is  slightly  different.  For  instance,  at 

8 .  =  85°,  it  results  in  a  null  near  the  specular  direction  for  the  smallest  roughness 

case.  When  the  analysis  was  extended  to  perturbation  method  regimes  (heights 

small  compared  to  a  wavelength)  the  behavior  is  quite  different.  In  that  case,  the 

narrowing  of  the  scattering  pattern  with  decreased  roughness  is  not  present.  For 

all  three  correlations,  the  magnitudes  decrease  with  decreasing  roughness  but  the 

individual  patterns  are  preserved.  The  Gaussian  and  power  law  patterns  are  broad 

(up  to  9  =  60°)  while  the  Bessel  function  correlation  shows  a  rapid  falloff  with 

increasing  scattering  angle  9  •  all  three  have  their  largest  values  for  small  8  .  To 

s  s 

summarize,  the  only  factor  that  introduced  any  variability  into  the  results  was  the 
use  of  the  Bessel  function  correlation  and,  even  for  that  case,  the  trends  are 
similar. 
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